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Abstract 

We show that the complexity of the Markov bases of multidimen- 
sional tables stabilizes eventually if a single table dimension is allowed 
to vary. In particular, if this table dimension is greater than a com- 
putable bound, the Markov bases consist of elements from Markov 
bases of smaller tables. We give an explicit formula for this bound 
in terms of Graver bases. We also compute these Markov and Graver 
complexities for all K x 2 x 2 x 2 tables. 

1 Introduction 

Let di,. . . ,dn be positive integers where di > 2. A multidimensional con- 
tingency table is an di X . . . X d„ array of nonnegative integers. Such a table 
represents the results of a census of individuals for which n discrete random 
variables Xi , . . . , Xn are observed (where we assume the random variable 
Xi takes values in [di] := {1, . . . ,di}) . Inferences about the collected data 
are made based on a statistical model or collection of models. This paper 
is concerned with the family of hierarchical log-linear models for which one 
assumes a set of interaction factors between the random variables [8J. Per- 
forming the exact test of conditional inference requires knowledge of the 
Markov basis of a given hierarchical model, which we describe below. 

When we assume that the sampling distribution of a table of observations 
is Poisson or multinomial, the sufficient statistics of any hierarchical model 
are given by certain marginal totals. The particular marginal totals that 
are sufficient statistics depend on the hierarchical model. For instance, for 
a di X ^2 X c?3 contingency table, the no three-way interaction model has 
sufficient statistics that are the three 2-way margins of the table: 

di d2 d-i 

i=l j=l k=l 
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where Uijk are the entries of the table. 

In general, a hierarchical model (and hence the marginal totals) is de- 
scribed by the list of the maximal faces Fi, . . . ,Fr of a simplicial complex A 
on n vertices. Computing marginal totals corresponds to a linear map from 
the space of tables to the space of marginals: 

r 

k=l 

where D = Yli=i di and = Wj(zp^ dj. This map is defined by 

r 

(^in,...,i„ : ij e [dj]) — > ^^ii,...,i„- 

Two tables t and u are said to be in the same fiber of vta if vta (t) = vta (u) . In 
other words, two tables are in the same fiber if they have the same margins 
with respect to A. We say that t and u are connected by the sequence 
of moves vi,...,Vs if each move Vj is in ker2(7rA) (that is, Vj has zero 
margins), t + X^iLi ^ table with nonnegative entries for each 1 < p < s, 
and u = t + Ya=i ^i- 

By a theorem of Diaconis and Sturmfels [H Theorem 3.1], for each hier- 
archical model given by A and d = (di, . . . , (i„), there exists a finite set of 
moves called a Markov basis such that any two tables that are in the same 
fiber of vta are connected by the moves in the Markov basis. Computing 
Markov bases via Grobner bases for the use in MCMC methods was initi- 
ated in [1], and since this first work computing Markov bases efficiently and 
describing Markov bases succinctly have been the major focuses of research. 
Recently substantial progress has been made. Simple Markov bases (consist- 
ing of moves with four nonzero entries) for decomposable models have been 
determined |5j, and similar Markov bases are known for reducible models 
[6l [9] . The case of binary graph models (where di = ■ ■ ■ = dn = 2 and A is 
a graph on n vertices) is worked out up to n = 5 [3]. 

The contribution of this paper is the most general form of a result first 
obtained in [1] for the no three-way interaction model for ii' x 3 x 3 tables. 
Our main theorem and its proof rely on ideas from [10] which treats the case 
of K d2 y- ■ ■ ■ y- dn tables where {2, 3, . . . , n} is a maximal face of A~ the 
so-called logit models. 

Theorem 1.1. Let A and d = {di , . . . ,dn) define a hierarchical model. Then 
there exists a constant m := m(A; d2, • • • , dn) such that for all di > m the 
universal Markov basis A4A,d consists of tables of the format r x ^2 x • • • x (i„ 
where r < m. 
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In other words, if we fix a hierarchical model together with the n — 1 
dimensions d2, ■ ■ ■ ,dn while varying the single dimension di , then for large 
enough di the universal Markov basis A4A,d will be obtained from Markov 
bases of small fixed-size tables; see Definition [3TTJ In particular, as a function 
of di, the complexity of computing and storing a Markov basis is bounded. 

We give the details of the proof of Theorem 11.11 in Sections 2 and 3, 
where we give an explicit computable upper bound for m{A; d2, ■ ■ ■ , dn)- 
We also present, in Section 3, a lower bound for m{A;d2,---,dn) which 
only applies in some cases. In the fourth section, we consider strengthenings 
and generalizations of the main result when more than one level is allowed 
to vary. In the final section we explicitly compute the complexity bound 
m(A; 2, 2, 2) for ah tables of the form K x 2 x 2 x 2. 

2 From Models to Matrices 

In this section we describe how to obtain a matrix A/\ corresponding to 
the linear transformation vta where the simplicial complex A describes a 
hierarchical model. We describe a decomposition for that is fundamental 
to the proof of Theorem 11.11 

Given the vector d = (di, . . . , dn), and a subset F = {ji < . . . < jg} C [n] 
we let dp = {djj^,dj^, . . . ,dj^). The columns of are in bijection with the 
D entries of a di x • • • x table, and we label each such column with the 
vector indexing the table entry (zi, . . . , in) € [di] x ■ ■ ■ x [dn]. Moreover, we 
order these columns lexicographically: 

(1,1,...,1,1) ^ (1,1,...,1,2) ^ ••• ^ (l,l,...,l,d„) -< 

(1,1,...,2,1) ^ (1,1,...,2,2) ^ ••• ^ (l,l,...,2,d„) ^ ••• 

(di,d2,...,d„_i,l) -< (di,d2,...,d„_i,2) ^ ••• ^ (di, (i2, . . . , 

Each row is labeled by a pair (F,e) where F = {ji,j2, ■ ■ ■ ,js} is a facet of 
A and e = (e^^, ej^, . . . , ej^) G [djj x [dj^] x ■ ■ ■ x [dj^] indexing the marginal 
corresponding to F. We first list the rows {F, e) where 1 S F. We impose a 
linear order on the facets where 1 G F and set {F, e) -< (G, /) if ei < /i, or 
if ei = /i and F ~< G, and in the case when ei = fi and F = G we use an 
arbitrary but fixed order of the indices. The rest of the rows will be listed 
again by some arbitrary but fixed order which will not play a role for the 
rest of the article. The entry of tIa in the column indexed by (ii, . . . ,in) 
and the row {F = {ji, . . . ,js}, {ej^ , . . . , e^J) will be equal to one if ij^ = Cj^ , 
ij^ = ej2, ■ ■ ■, and ij^ = Cj^; and it will be zero otherwise. 



3 



Example 2.1. Let A = {{1,2}, {1,4}, {2, 3}, {3, 4}} and di 
di = 2. This is the binary 4-cycle model. 



d2 = ds 
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Here the 16 columns are indexed as (1, 1, 1, 1) -< (1, 1, 1,2) -< (1, 1, 2, 1) -< 
(1,1,2,2) -< ••• -< (2,2,2,1) -< (2,2,2,2). The first four rows arc indexed 
by (Fi, (1, 1)), (Fi, (1, 2)), (F2, (1, 1)), and (F^, (1, 2)) where Fi = {1, 2} and 
F2 = {1,4}. The second block offour rows are indexed by (Fi, (2,1)), (Fi, (2,2)), (F2, (2,1)), 
and (F2, (2, 2)). For the rest of the rows we have chosen the order (F, (i, j)) -< 
(G, {s,t)) if F = {2,3} and G = {3,4}, or if F = G and ~< {s,t) lexi- 
cographically. 

We observe that when the rows and columns are ordered as described, the 
matrix A/\ exhibits a block structure. In the above example, the upper-left 
and lower-right blocks of the first eight rows are identical augmented with 
the two blocks of zeros. And the last eight rows are split into two identical 
matrices. We summarize this observation in the following lemma where we 
assume the ordering of the columns and rows of that we introduced 
above. 

Lemma 2.2. Let A = {Fi, . . . , F^} and d = (di, . . . , d„) define a hierarchi- 
cal model. Then 



A 
0^0 



B B B 






A 
B 



where A is a Ylk=ii-^k/di) x (D/di) matrix with Fi 
containing the vertex 1, and where B is a 
Hence there are di copies of A and B . 



s+l 



■ . . , Fg being the facets 
Dfe X (D/di) matrix. 
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Remark. The matrices A and B are also matrices that come from hierarchical 
models. Note that the matrix A is the matrix A^ for the simplicial complex 

r = link(A) := link(A, {1}) = {F \ {1} | F G A and 1 G F} 

and the vector d' = (2, 3, . . . , n). The matrix B is the matrix for the 

simplicial complex 

A \ {1} = {F I F G A, 1 ^ F, and F is a facet}. 

3 Proof of the Finiteness Theorem 

In this section, we proceed with the proof of Theorem II. 1[ To do this, we 
will prove a finiteness theorem for the Markov bases of arbitrary matrices 
which come in a block form akin to the one demonstrated in Lemma 12.21 

Definition 3.1. Let A G N'^^"' be an integer matrix with no zero columns. 
A finite set A4 C kerz(A) of integer vectors in the kernel of A is called a 
Markov basis of A if any two nonnegative integer vectors in the same fiber 
of A can be connected by a collection of the elements in Ai. That is, for 
any t, u G with At = Au, there exists a sequence of moves {vj}|^-^ C M 
such that 

p s 
t + Vj > for all 1 < p < s and t + Vj = u. 

i=l 1=1 

A Markov basis of j4 is called minimal if no subset of is a Markov basis 
of A. The universal Markov basis A4{A) of A is the union of all minimal 
Markov bases of A. 

When A = Aj\ is the matrix associated to a hierarchical model, we use 
the shorthand MA,d to denote the universal Markov basis of Aj\. 

Definition 3.2. Let u,v, and v' be nonzero vectors in kerz{A). We say 
u = V + v' is a conformal decomposition of u if > implies < Vi,v^ < Ui, 
and Ui < implies Ui <Vi,v[ < for all 1 < i < n. The set Q{A) C ker2(A) 
of integer vectors with no conformal decompositions is called the Graver 
basis of A. 

One can show that Q{A) is a finite set [111 Chapter 4] and any minimal 
Markov basis of ^ is a subset of G{A) [TT, Chapter 5]. Thus, M{A) is a 
finite set. 

Definition 3.3. Let ^4 be a d x n matrix with columns oi, . . . , a„ and B 

be a p X n matrix with columns 6i, . . . , 6„. The r-th generalized Lawrence 
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lifting of A with B is the {rd + p) x rn matrix A(A, B ,r), whose columns 
are the vectors 

A{A, B, r) = {ai(S> Cj (Bbi\l <i < n,l < j < r}. 

In particular, the matrices A/\ are of the form A{A, B, r) with A = 
^iink{A)) B = j4^\^{i} and r = di. When B is the n x n identity matrix and 
r = 2, the matrix A{A, B, r) is called the Lawrence lifting of A [11, Chapter 
7]. For B = In but general r, this matrix is called the rth Lawrence lifting 
of ^ [TU]. 

Remark. An integer vector in the kernel of A(^, B, r) can be represented as 
an r X n matrix where each row is in the kernel of A, and the sum of the 
rows is in the kernel of B. For instance, the following 2x8 matrix is the 
representation of such a vector in the kernel of j4a in Example 12.11 

"2-2-1 1-2 2 1 -1 " 
-1 1 1-1 1-1-1 1 ■ 

Definition 3.4. The type of a vector in represented as an r x n matrix is 
the number of nonzero rows of this matrix. The Markov complexity m(A, B) 
of a d X n matrix A and apx n matrix B is the largest type of any vector in 
the universal Markov basis of A{A,B,r) as r varies. Similarly, the Graver 
complexity g{A, B) of these two matrices is defined as the largest type of 
any Graver basis element of A{A,B,r) as r varies. Analogously, we define 
g{A;d2,---,dn) and m{A;d2,---,dn), the Graver and Markov complexities 
of the hierarchical models corresponding to A as di varies. 

We will show that the Graver complexity g{A, B) is finite. This implies 
that the Markov complexity is also finite since m{A, B) < g{A, B). In order 
to do this we relate the Graver basis G{B ■ G{A)) to the collection of Graver 
bases Q{A{A, B,r)). We emphasize the "double" Graver construction: we 
first compute the Graver basis of A and obtain the set GiA). We consider 
each element in G{A) as a column vector. Then the vectors B ■ G{A) are 
computed by multiplying each element of G{A) with B. Thus, B ■ G{A) is a 
p X 1^(^)1 matrix. Finally we compute the Graver basis of B ■ G{A). 

Theorem 3.5. The Graver complexity g{A,B) is the maximum 1-norm of 
any element in the Graver basis G{B ■ G{A)). 

In order to prove the above theorem we need the following lemma. 

Lemma 3.6. Let u = [u^;u'^; . . . be in the Graver basis of A(A,B,r). 
Suppose that = u ^ + v"^ is a conformal decomposition where and v"^ are 
in the kernel of A. Then the element [u^; . . u^~^; u^; v^; n*"*"^; . . . ; u^] is in 
the Graver basis of A{A, B,r -\- 1). 



6 



Proof. Suppose not. Then [u^ ; . . . ; u"^ ^; ; v"^ ; u^^^ ; . . . ; li^] has a conformal 
decomposition 

\u^;...]v^~^;v^]v'^;ii^+^]...]v'] + [u^ ; . . . ; u'-^ ; ; ; u'+^ ; . . . ; u'] 

where both vectors are in the kernel of A(A,B,r + 1). Now since = 
+ v'^ = {v^ + v^) + (v^ + v'^) is a conformal decomposition of n', so is 
(v^ + v^) + {v^ + u^). We note that neither the first nor the second sum is 
zero. But then the two nonzero vectors [u^; . . . ; u^^^; + -u^; u^^^; . . . ; n*"] 
and [n^; . . . ; u^~^;v^ + u^; . . . ; u"^] are in the kernel of A{A, B, r), and 
their sum forms a conformal decomposition of u. This is a contradiction 
since u is in the Graver basis of A(A, B, r). □ 

Proof of Theorem \3.5\ Lemma 13.61 implies that in order to compute the 
Graver complexity g{A, B) we only need to consider elements u = [u^; . . u^] 
where is in the Graver basis G{A) = {vi, . . . Given any such u, we 

construct a vector T G Z'^ where the ith entry counts how many times Vi 
appears in u. The 1-norm of T is the type of u. Hence we need to show 
that r is in the Graver basis of B ■ G{A) if and only if u is in the Graver 
basis of A{A,B,r). If T is not in the Graver basis of B ■ G{A), then it has 
a conformal decomposition Fi + r2 such that B ■ GiA) ■ Tj = for i = 1,2. 
Reversing the operation that produced F from u, Fi and F2 yield vectors 
vi, V2 £ ker2(F(A, i?, r) such that u = vi + V2 and this decomposition is 
conformal. Thus u could not be in the Graver basis of A(A, B ,r). 

Conversely, a conformal decomposition of u translates into a conformal 
decomposition of F since none oi u^, . . . G G{A) have conformal decom- 
positions. □ 

Proof of Theorem The hierarchical model defined by A and d gives rise 
to A/\ which is of the form A{Ai^^^^-j, A^\^^i^, di) by Lemma |2.2[ Theorem 
13.51 implies that the Markov complexity ?Ti(^iink(A); ^A\{i}) is bounded by 
the finite Graver complexity 5'(^iink(A)7 ^A\{i})- This means that for all 
di > "i'(Aink(A)) ^A\{i}) the universal Markov basis M.A,d will consist of 
tables of the format r x ^2 x • • • x d„ where r < m{Aliy^^^^, A^\^iy). □ 

In practice, the Graver complexity and the Markov complexity may vary 
a lot, as the following examples illustrate. All of our examples were com- 
puted using 4ti2 [7J and the results for Markov bases of reducible models 
using P |9| . 

Example 3.7. Let A = {{1, 2}, {1, 3}, {2, 3}}. Then for ^2 = ^3 = 3, 
the Markov complexity is m(A;3,3) = 5 (the main result in [1]), while the 
Graver complexity is 5(A;3,3) = 9. 
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1 2 




3 3 

Figure 1: Renaming the vertices in A can change the Markov complexity 

Example 3.8. Let A = {{1, 2}, {1, 3}, {1, 4}, {2, 3}, {3, 4}}. Then for ^2 = 

d^ = d4 = 2 the Markov complexity is m(A; 2, 2, 2) = 2 while ^(A; 2, 2, 2) = 
4. On the other hand, for the same complex in a different orientation (see 
Figure 1) A' = {{1, 2}, {1, 3}, {2, 3}, {2, 4}, {3, 4}} and ^2 = ^3 = ^4 = 2, 
the Markov complexity is m(A'; 2, 2, 2) = 4 while ^(A'; 2, 2, 2) = 16. 

The Graver complexity g{A, B) gives an upper bound for the Markov 
complexity m{A,B), in terms of the Graver basis of B times the Graver 

basis of A. There is an analogous lower bound for the Markov complexity 
in terms of the Graver basis of B times the Markov basis of A. To describe 
this lower bound, we introduce the notion of semiconformal decompositions. 

Definition 3.9. Let u, v, and v' be nonzero vectors in \cti{A). Wc say that 
u = V + v' is a semi-conformal decomposition if > implies that Vi < Ui 
and v'i < implies Ui < v'i for all 1 < i < n. Note that if the first condition 
holds {vi > implies that Vi < Ui for all i) the second condition { v'i < 
implies Ui < v'i for all i) is satisfied automatically. The set S{A) C keri{A) 
is the set of vectors in kerz{A) which have no semi-conformal decomposition. 

A useful fact about vectors in the kernel of a matrix which have no semi- 
conformal decomposition is that they must belong to the Markov basis. 

Lemma 3.10. Let A G I^'*^" with no zero columns. If M. is any Markov 
basis of A, then S{A) C M. In particular, S{A) C M{A). 

Proof. Suppose that u € S{A) has no semiconformal decomposition but 
there is some Markov basis Ai oi A that does not contain u. Write u = 
u"^ — u" as the difference of two nonnegative integer vectors with disjoint 
support. Note that u"*" and u~ belong to the same fiber. Since is a 
Markov basis for A there is a sequence of elements from M, {v^, v^, . . . , v*} 
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with s > 2 which connects u to u+ where intermediate summands are 
always nonnegative. In other words, we can write 



and the set of indices where v is negative is a subset of the set of indices 
where u~ is nonzero, and u~ > \vl\ for this subset of indices. But this 
imphes that 



k=2 

is a semiconformal decomposition of u. This contradicts our assumption 



Theorem 3.11. The Markov complexity m{A,B) is bounded below by the 
maximum 1-norm of any element in the Graver basis Q{B ■ S{A)). 

Proof. Let T € Q{B ■ S{A)). FoUowing the proof of Theorem 13.51 T trans- 
lates into a vector u = [n^; . . . with each ti* € S{A). Furthermore, we 
know that u lies in the Graver basis of h.{A,B,r). We wish to show that it 
lies in some minimal Markov basis of A(A, B, r). To do this, we show that u 
has no semiconformal decompositions. Suppose, to the contrary that there 
was some semiconformal decomposition of u. Since u is in the Graver basis 
of A(yl, i?, r), any semiconformal decomposition u = v + v' induces a semi- 
conformal decomposition of (at least) one of the vectors n*. However, this 
is a contradiction, because 5(^4) consists of vectors with no semiconformal 
decompositions. Thus r, which is the 1-norm of F is a lower bound for the 
Markov complexity m(A, i?). □ 

Example 3.12. Applying Theorem 13.111 for the simplicial complex A = 
{{1, 2}, {1, 3}, {2, 3}} we have m(A; 3, 3) > 5, m(A; 3, 4) > 8, and m(A; 3, 5) > 
12. It is interesting to note that in the first two of these cases our lower 
bound for the Markov complexity is equal to the value reported in [U [2]. 
The Markov complexity M(A; 3, 5) remains undetermined. 

4 Generalizations and Extensions 

Given Theorem II. H it is natural to ask to what extent this result is the 
best possible. In particular, do there exist any bounds on the complexity of 
Markov basis elements if we fix A and allow di and d2 to vary? The answer 
to this question is negative if we allow arbitrary A. 



s 




s 




that u G S{A) and hence u G A4. 



□ 
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Example 4.1. Let A = {{1, 2}, {1, 3}, {2, 3}} and fix ds > 2. There are no 
bounds nil and m2 such that every Markov basis element has format strictly 
contained in a mi x m2 x ^3 table. A well-known example of such a move 
of large format is shown in ^y. Denote by eijk the 3- way table with a 1 in 
the {i,j,k) position and zeroes elsewhere. For each m> 1 the vector 

U = em + 6221 H h emml 

+6122 + 6232 H h 

— 6112 — 6222 — • • • — 6mm2 

— 6121 — 6231 — • • • — Cm-lml — Cmll 

belongs to the universal Markov basis M{Aj\) that has format m x m x 2. 

On the other hand, for reducible models, we can prove more general 
finiteness results. These are based on the structural theorems for building 
Markov bases for reducible models in [6l [9]. Recall that for a simplicial 
complex A, |A[ = UpeAF is the underlying set of A. 

Definition 4.2. A simplicial complex A is called reducible with decompo- 
sition (Ai,5, A2), if 5 € Ai, 5 G A2, |Ai| n IA2I = S, and A = Ai U A2. 

If A is reducible we use the notation d^ and to denote the substrings 
of d that are indexed by |Ai| and IA2I respectively. 

Example 4.3. The simplicial complex A = {{1, 2}, {1, 3}, {2, 3}, {2, 4}, {3, 4}} 
pictured in Figure 1 is reducible with S = {2, 3}, Ai = {{1, 2}, {1, 3}, {2, 3}}, 
and A2 = {{2, 3}, {2, 4}, {3, 4}}. The vectors d^ and d'^ are (^1,^2,^3) and 
{d2,d3,d4), respectively. 

One of the main results from [6l [9] is the constructive version of the 
following. 

Theorem 4.4. Let A be a reducible simplicial complex and let d be given. 
Letli andl2 be the maximum 1-norm of any element ofM.^^,11 and J^i^^,p, 
respectively. Then the maximum 1-norm of any element of A4.A,d is max{4, /i, I2}. 

This allows us to deduce that reducible models have Markov bases of 
finite complexity as many levels vary. 

Corollary 4.5. Let A be a reducible simplicial complex with induced sub- 
complexes and A2 and suppose that 1 € |Ai|\|A2| and 2^ |A2|\|Aij. Let 
ds, . . . ,dn be given. Then there exists constants (mi, m2) = m(A; ^3, . . . , d„) 
such that every element in the universal Markov basis TWA.d has format 
smaller than mi x m2 x d^ x ■ ■ ■ x dn- 
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Proof. Restricting to Ai and A2 and allowing di and ^2 to vary respec- 
tively, we know by Theorem 11.11 there is a bound on the format of Markov 
basis elements that appear in M/^-^ and d?- However, a bound on 
the format also implies that these vectors have bounded 1-norm (one such 
bound is the Markov complexity times the largest 1-norm of any element in 
^(^iink(A)))- Applying Theorem 14.41 we deduce that every element of A^A,d 
has 1-norm bounded by some fixed constant. But bounded 1-norm implies 
bounded format and completes the proof. □ 

Besides the condition that A is reducible, the crucial requirement to 
prove the preceding corollary was that 1 and 2 were not adjacent to each 
other in A. We conjecture that this property is enough to guarantee bounded 
Markov complexity in general. 

Conjecture 4.6. Suppose that {1, 2, . . . , j} is an independent subset of the 
underlying graph of A. Then for fixed dj^i, . . . ,dn, there exists numbers 
(mi, . . . , ruj) = m(A; dj+i, . . . , dn) such that every element in the universal 
Markov basis MA,d has format smaller than mi x • • • x mj x dj+i x ■ ■ ■ x dn. 

Unfortunately, we do not know if Conjecture 14.61 is true even in the sim- 
plest nonreducible case, namely the four-cycle A = {{1, 3}, {2, 3}, {2, 4}, {1, 4}} 
with j = 2. If a proof exists, it must depend on techniques different from 
those developed here, because if j > 1 and A satisfies the hypotheses of Con- 
jecture 14.61 there is no bound on the formats of the Graver basis elements 
for ^A- 

5 Computations 

The following table displays computational results of the Markov complex- 
ity and Graver complexity of all binary hierarchical models where one of 
the dimensions of the tables is allowed to vary. In the notation of Theorem 
11.11 this is the Markov complexity m(A; 2, 2, 2) and the Graver complexity 
(/(A; 2, 2, 2). Note that all the entries which are marked with a star are 
Markov and Graver complexities which were not known, or could not have 
been determined, without the use of Theorem 11.11 All of the computations 
described in this section were performed using the toric Grobner basis pro- 
gram 4ti2 [7|. The second and fifth column m correspond to the computed 
Markov complexity and the third and sixth column g is the Graver complex- 
ity. We use the bracket notation from multivariate statistics [8] for denoting 
simplicial complexes. Thus [12] [23] [34] represents the simplicial complex 
{{1,2},{2,3},{3,4}}. 
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Model 


m 


5 


Model 


m 


5 


[123] [124] [134] [234] 


2 


2 


[12] [13] [23] [34] 


2 


10* 


[123] [124] [134] 


2 


2 


[12] [13] [23] [14] 


2 


2* 


[123] [124] [234] 


2 


2 


[12] [23] [24] [34] 


4 


8* 


[123] [124] [34] 


2 


2 


[12] [14] [23] [34] 


4* 


5* 


[123] [234] [14] 


4 


4 


[123] [4] 


2 


2 


[123] [14] [24] [34] 


4* 


4* 


[234] [1] 


2 


8 


[234] [12] [13] [14] 




12* 


[12] [13] [23] [4] 


2 


8* 


[12] [13] [14] [23] [24] [34] 




10* 


[23] [24] [34] [1] 


4 


16* 


[123] [124] 


2 


2 


[12] [23] [34] 


2 


6* 


[123] [234] 


2 


2 


[12] [14] [23] 


2 


4* 


[123] [24] [34] 


4 


4* 


[12] [23] [4] 


2 


6* 


[234] [12] [13] 


2 


10* 


[12] [13] [4] 


2 


3* 


[123] [14] [24] 


2 


2* 


[23][34][1] 


2 


14* 


[12] [13] [23] [24] [34] 


4 


16* 


[12] [34] 


2 


4 


[12] [13] [14] [23] [24] 


2 


4* 


[12] [3] [4] 


2 


4* 


[123] [34] 


2 


2 


[34] [1] [2] 


2 


12* 


[123] [14] 


2 


2 


[1][2][3][4] 


2 


10* 


[234] [12] 


2 


4 









Finally, we used the Theorem 13.111 to compute some lower bounds for 
the Markov complexity. 

m([12][13][23],3,5) > 12 

m([12][13][23],4,4) > 16 

m([123][124][134][234],3,3,3) > 19 

These lower bounds are benchmark values for extending the types of 
results pursued in [T| and [2] in which the values m([12] [13] [23], 3, 3) = 5 
and m([12][13][23], 3, 4) = 8 were explicitly computed. 
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